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ABSTRACT 


Knowledge of the elastic properties of composite materials can be an 
invaluable tool for both the quality assurance of manufacturing techniques and 
design verification. Recent advancements in ultrasonic velocity measurements 
have demonstrated the ability to recover elastic properties in anisotropic 
laminates. A simplified experimental setup was investigated to recover the 
elastic properties based upon the flexural wave propagation in anisotropic 
laminates. The initial objective of this thesis was to verify flexural wave 
propagation in composite laminates through the comparison of experimental and 
theoretical phase velocities. In the second part of this thesis, the experimental 
phase velocities were used to calculate the elastic properties of the material by 
inverting the governing equations. The initial method used to recover elastic 
constants was successful in the recovery of a partial set of the bending and 
extensional stiffnesses. The inability to recover all bending stiffnesses dictated 
the investigation of a second method. This method used an iterative method 
based upon a nonlinear Newton's method to recover the bending stiffnesses. This 
method did not converge due to ill conditioning of the solution matrix. Although 
this method did not converge, it is believed that other more robust methods 


suggested herein would converge to the proper solution. 


lil 


> Y 
i (— 


y 7 Fe 
Aa: 


TABLE OF CONTENTS 


I. INTRODUCTION Ph 'icccicccccccccage Sorte e cere fo a seen l. 

[Il WAVE THEORY FOR THIN®PLATES i .c2cce cee ee 3. 
A. CLASSICAL PLATE THEORY sein ..<......-0¢2:-c eee ee De 
B. HIGHER ORDER PLATE THEOR Y (45"...... cogureees en 8. 

Il]. DISPERSION MEASURE MIEBIN TS rr eeter eee seen eee 14 


B. EXPERIMENTAL RESULTS... este 16 
IV. ALGORITHMS FOR RECOVERING MATERIAL CONSTANTS......... oF 
A. (O}ig LAMINATE... serene cent D7. 

B. [0,45,90)-45lo5 LAMINATE. 32) sees re BD. 

C. [0,90l46°LAMINATE........225..,. ee ee one 

D. SUMMARY/ DISCUSION OF ALGORITHMS USED FOR RECOVERY 

OF ELASTIC’ CONSTANTS....00)........enen 42. 
V. SUMMARY... oo... .cc..cccssessesthes ss) een oe ee 45. 
APPENDIX A- LAMINATE MATERIAL PROPERTIES........0cccccscsceee 47. 
APPENDIX B- EXPERIMENTAL DATA =2).....::9 0s ee eee 49 
APPENDIX C- EXPERIMENTAL ERRORS .Q.u...cccccccccscssscessscssssessecsseeesees 56. 
LIST OF REFERENCES. 225 ee 59. 
[INITIAL DISTRIBUTION LIST. ...............0 ee 61 


I. INTRODUCTION 


Composite materials have been used quite extensively in the past for high 
performance structural components and promise to be major structural materials 
of the future [Refs. 1-4]. The inherent advantage of composite materials is their 
high strength to weight ratio. Composite structures present unique problems in 
that these materials are produced by complex manufacturing processes. These 
processes include a large number of manufacturing variables which increase the 
number of opportunities for introduction of defects. Variables such as state of 
cure of the resin, fiber to matrix bond, and fiber to resin ratio have an effect on 
the elastic properties of a laminate. It is thus imperative that controlled 
procedures be established for quality assurance throughout the manufacturing 
process. A knowledge of the elastic properties of composite materials would be 
invaluable for both design verification and quality assurance. There is a need for 
a simple and quick in situ test which can yield quantitative parameters descriptive 
of the mechanical state of the material. 

Several non-destructive testing (NDT) techniques with their own inherent 
strengths and weaknesses are currently available to assess the mechanical states of 
composites and are discussed below. 

Radiographical techniques have long been available to detect large voids, 
porosity content, foreign debris and other imperfections in composite laminates. 
The presence of these imperfections has a most definite impact upon the material 
properties. There seems to be disagreement on the ability of radiography to 
determine fiber volume while the value of this technique reduces with increased 


laminate thickness. [Ref. 4] 


Ultrasonic methods are one of the most widely employed techniques for the 
non-destructive testing of composite laminates [Ref. 5-10]. Ultrasonic NDT 
techniques include velocity measurement, ultrasonic spectroscopy and 
attenuation measurement. 

There is a definite relationship between void content and mechanical 
strength. Ultrasonic spectroscopy has proven itself to be a viable technique to 
detect and locate voids, delaminations and other defects. Also, the dependence of 
ultrasonic attenuation upon void content has been clearly demonstrated. While 
these NDT techniques are effective tools for testing of production quality 
control, they can not quantify the actual material mechanical properties. [Ref. 4] 

Several works on ultrasonic velocity measurement to recover the material 
properties have been published recently. Dayal and Kinra [Ref. 5] verified the 
Lamb dispersion for an anisotropic laminate immersed in a water bath by 
measuring the phase velocities and attenuation of ultrasonic waves in thin plates. 
Every and Sachse [ Ref. 6] have shown through numerical simulation it is 
possible to recover a partial set of elastic constants of a mildly anisotropic 
material using the quasi-longitudinal mode. However, the quasi-longitudinal 
mode is typically difficult to excite in thin plates. Wu and Ho [Ref. 7] measured 
the energy (group) velocity of the quasi-longitudinal and quasi-shear mode 
through the thickness of a 36 ply unidirectional laminate using an ultrasonic 
technique. The energy velocity measurements were used to calculate the phase 
velocities through a numerical method. From the phase velocities the elastic 
constants could be calculated. Wu and Chiu [Ref. 8] found good agreement 
between theoretical and experimental results when exciting the shear mode in a 


unidirectional laminate. When attempting to determine the elastic constants from 


tO 


velocity measurements of the shear mode, they reported that even a small 
perturbation in the shear wave velocity will induce large differences in the 
calculated elastic constants. Castagnede et. al. [Ref. 9] utilized the nondispersive 
bulk modes to determine the elastic constants of a glass/ epoxy unidirectional 
composite. Every and Sachse [Ref. 10] used the group velocities of the bulk 
longitudinal mode to recover the elastic constants of an anisotropic solid. Every 
and Sachse utilized the point-source--point-receiver (PS-PR) technique to 
measure the group velocities. The group velocities were used to determine the 
elastic constants of an anisotropic solid. The PS-PR technique requires access to 
both sides of the specimen which is not always possible when measuring in situ. 

While these methods have shown success in recovering material properties of 
composite laminates, there are limitations upon their usefulness to the composite 
community. Castagnede et. al. and Every and Sachse have proven successful in 
recovering elastic properties while using the bulk modes excited by the PS-PR 
technique. However, excitation of the bulk mode is difficult in thin plates. Wu 
and Ho and Wu and Chiu successfully utilized immersional ultrasonic 
measurements to recover the elastic constants. One drawback of this method is 
that the test specimen must be immersed in water. 

This thesis investigates the possibility of using phase velocity measurements 
based on plate wave modes in thin composite laminates. The method utilizes a 
single pulser and receiver combination on the same side of the plate. The 
experimental setup does not require immersion of the specimen nor does it 
require access to both sides of laminate. The experimental setup is discussed in 


further detail in Chapter III. 


Gorman [Ref. 11] has shown for this type of transmitter/receiver setup, 
when the wavelength is much larger than the plate thickness, only two modes of 
wave propagation are present with a significant amount of energy, the lowest 
order extensional and the flexural modes. Initially, this thesis is concerned with 
the verification of the dispersion curves predicted by higher order plate theory 
of [Ref. 12]. Experimentally measured phase velocities are compared with 
theoretical phase velocities calculated using the manufacturer's supplied material 
properties. Next, the possibility of recovering the material properties is 
explored through the utilization of both the classical plate theory and higher 
order plate theory. I[t is shown that by inverting the governing equations for the 
flexural mode, the experimentally measured phase velocities can be used to 


recover a partial set of the material properties of several experimental laminates. 


Il. WAVE THEORY FOR THIN PLATES 


In the experimental section to follow, the experimentally measured phase 
velocities are compared to theoretical predictions. Therefore, the theory of 
wave propagation in thin plates must be understood. In this section, the classical 
plate theory and the higher order plate theory of wave propagation are 


presented. 


A. CLASSICAL PLATE THEORY 
To understand the propagation of flexural waves in thin plates, we begin by 
considering a plate of thickness h and of infinite extent in the x-y plane as 


depicted tn Figure 2.1. 





Figure 2.1 Plate Element (after Graff [Ref. 13]). 


A differential element h dxdy will be subject to the various shear forces, 
bending and twisting moments as shown due to the applied stresses. Summing 


forces in the z - direction gives 


IQx IQy O°w 
=O) ; 72, 1\ 
a Oy es ore 2-1) 











W(x, y,t,) 18 the displacement of the mid-plane of the plate in the z-direction. 


Similarly, summing moments gives 











OM, OMy 
. O70 NP 
Oy OX ; og 
a ae (2.3) 
OX Oy 


Solving (2.2) and (2.3) for the shear forces Q, and Qy and substituting into 


(2.1) results in a single governing equation in terms of the various moments, 








2 2 2 2 
oh a 2 ee (2.4) 
Ox“  dxdy odydx dy ot 


Recognizing Mxy = -Myx, (2.4) simplifies to 











OM, OM OM O*w 
2 ==) Fi : 2.5 
ax?” Oxoy i oy* se (2.9) 














For a laminate arranged symmetrically about its mid plane, the bending 


moments are related to the curvatures by 


Mx Di: Di2 Die || Kx 
My |=|Di2 De22 Dee] Ky |. (2) 
Muy D6 D 26 Des Kxy 


The bij are the bending stiffnesses and KxKy and Ky are the middle plane 


curvatures which are defined as 





Kx a 
W 
Ky |= oye : (2) 
Kxy O2w 
 axay | 


The displacement of the flexural wave is given by 
a= Acwiuga cae!) (2.8) 

The value k is the wave number, |; and lg are the direction cosines in the x 
and y directions, respectively, @ 1s the circular frequency and A is the amplitude 
of the wave. 

Substituting (2.6) through (2.8) into (2.5) and setting the body forces to zero 
gives 

k*(h* Du + 4h beDie + 2h? (Diz + 2Do6) + 4hil2’ Dae + 12* D22) = ph’. (2.9) 

For flexural waves propagating in the x direction (i.e. lj =1 and l9=0), 

equation (2.9) reduces to 
k*D11 = pha’. (2.10) 

Similarly, for flexural waves propagating in the y direction equation (2.9) 

reduces to 
k*D22 = pho’. lL) 
From elementary wave theory, the phase velocity, Vp, in terms of circular 


frequency and wave number can be written as 


Ww 
V, = te jag ee 
re (2.12) 
Substituting (2.12) into (2.10) and solving for the phase velocity of waves 


propagating in the x-direction gives 


21 
Ve calle )4 = (oceee 


a7 6: 
Likewise, substituting (2.12) into (2.11) and solving for the phase velocity of 
waves propagating in the y-direction yields 
D22w* 


. 4 r2Do2f2 - 
= : 2.14 
p= (SEO yt = (Aree (2.14) 


It can be observed that the phase velocity is a function of the frequency, Le., 
the phase velocity is dispersive. For low frequency conditions classical plate 
theory predicts the correct response. However, for short wavelength, high 
frequency conditions unbounded phase velocities are predicted. This physically 
unreasonable condition is the result of neglecting rotary inertia and shear effects 


in the theory. 


B. HIGHER ORDER PLATE THEORY. 

Classical plate theory is based on the Kirchoff hypothesis whereby the 
transverse shear deformation effects are neglected. However, the shear 
deformation effects are significant for a composite laminate owing to the fact 
that a composite laminate has a low transverse shear modulus. Consequently, the 
effects of transverse shear must be included to accurately predict the physical 
condition. The dispersion theory presented by Tang, Stiffler and Henneke [Ref. 
12] for composite plates is an extension of the work of Mindlin [Ref. 14] 
performed on isotropic plates, which accounts for both transverse shear and 
rotary inertia Corrections. 

We begin our discussion of higher order plate theory by considering a 
composite laminate plate of thickness h with the coordinate system as presented 


in Fig. 2.1. The following displacement fields are assumed: 


U=Uo(X,V.t)+ Zwx(X, yt) (2215) 


V =V0o(x, Y,t)+Zw(Xx, y,f) 26) 


W=W(X,y,L). CAG 

The displacement coordinates u, v, and w are in the x, y, and z directions, uo 
and Vo are the mid-plane displacement components, yxand wy are the rotation 
components along the x and y axes, respectively. Substituting (2.15) through 


(2.17) into the small strain-displacement relations of elasticity theory yields : 














OUo OW« 5 
poe 2h 
E 5 +Z 5y (als) 
OVo OWy 
ae 22S, 
y ay ae ay ( ) 
E2=0 2.20) 














The stress strain relations for a orthotropic lamina, including the out of 


plane shear terms, are 
Ox Ore Ore 0 0 Qie || Ex 


Oy Qi2 Qe 0 O Qazoell Ey 
=| Q 0 Qa Qs O |] Yi. ee) 

t| |0 0 Qs Qs 0 | % 

Try | | Qe Qe 0 0 Qes | Voy 
Qij for i=1,2,6 are plane stress reduced stiffnesses, and Qij for 1=4,5 are 


transverse shear stiffnesses. The force and moment resultants per unit length 


acting on the laminate are obtained by the integration of the stresses in each 


lamina through the plate thickness, 


hi2 
(Nx,Ny,Ny)= | (Ox,Oy, Tz O25) 


—h/le2 
hl2 


(Qx,Qy) = | (Te, Gea? (2.26) 


-h/2 
h/2 


(Mx,My,Mxy) = | (Ox,Oy, Ty»zaz Og5)) 


~hi2 


Substituting equations (2.18-23) and (2.24) into equations (2.25-27) gives 





QUo 

Ox 

OVo 
Nx Au Az 0 O Aw Bu Bie Bie dy 
Ny Ai2 Azz OQ O Az Bi2 Bo Bae CME hy, 
Q||0 0 Aw As 0 0 0 off % 








My Bi2 Boe Bos Di2 D2 Dee oon Cas 
May Bis Bee Bes Die Dae Dee OWy 
dy 
ay Ox 


The extensional stiffness, Aij, coupling stiffness, Bij, and bending stiffness, 


Dij, are given by 


hl2 


(Aj,By,Dy) = | (Qy)k(1,2,2°)dz ——i,j=1,2,6 (2.29) 
-h/2 
and 
hi2 
Ai = KK; | (Qy)xdz i,j=4,5 (2.30) 


-h/2 
The shear correction factors are included to account for the fact that the 


transverse shear distributions are not uniform across the thickness of the plate. 
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Summing forces and moments on a differential element and neglecting the 


body forces, the equations of motion are reduced to: 









































ON« ONxy . 0" Uo O° Wx 
—— + = +A re | 
Ox oy Por ot? aie 
2 2 
ON xy a ONy _ 0 O = O ua (se) 
ax oy at at 
O Ox OQny a O =e (2.33) 
Ox Oy ot 
OMx OM "Uo O Wx 
= Qx = R / 4 
Ox a ot® ot® 2.58) 
OMny OM, O°Vo O Wy 
-Q,=R / 23) 
axa ge AE 2.39) 


P is the mass density and 


hi2 
(pA)= | p(t.z,z°)az. 


-hi2 


(2.36) 


By substituting equation (2.28) into equations (2.31-35), the equations of 
motion in terms of the displacements and rotation are obtained. 

For this thesis, only symmetric laminates are considered. This results in 
considerable simplification because the coupling stiffnesses, Bij, and the normal- 
The equations of 


rotary inertia coupling coefficient ,R, are identically zero. 


motion for the flexural mode reduce to 
2 2 2 2 
OWx O°W Fog is MDD pA (ipa 5 “2 (2.37) 


























A + ap ; 
2 
peep.) +n 2 popes 20. 
Ds Oxoy ay Ox OxOy 
2 
Saye = )o ; (2.38) 


1] 




















O° Wx OW “ ony o* Ow 
DiS + 2D + De 5 + (Die Dos) — 26 5/2 Asle tS) 
OW O° wx 
= —w- —/ =) 2.39 
al la — 


For flexural wave propagation , we consider plane waves described by 


wis VER Sols, (2.40) 
Wx Be Tp (2.41) 
Wy as yee oy (2.42) 


The value k is the wave number, |! and 12 are the direction cosines of the 
wave vector in the x and y directions, @ is the circular frequency, and W, yx 
and yy are the amplitudes of the plane waves. After substituting equations 
(2.40- 42) into equations (2.37 - 39), the determinant of the resulting set of 
equations gives the dispersion (characteristic) equation for flexural wave 
propagation. Because this study is limited to symmetric quasi-isotropic 
laminates, in addition to Bij and R=0, Aig = Avo = Ags = 0 and Dig =D 26. This 


results in the following characteristic equation for flexural wave propagation. 


DiikK*le +2Diek*hl2  Diek® +(Di2+ Des) K 7h 2 lAsskl; 
+DesK*l2° + Ass—lw* 
Disk*® +(Di2+Des)k*hle = Desk*li? + 2D16K “Il 2 iAaaklo =0 
+D22k*lo+ Ass —lw* 
| —iAsskl; ~iAsskle —Assk?1? — Aask?I2? — p' a? | 


(2.43) 
Several simplifications can be made for waves propagating along principal 
directions. For waves propagating in the x- direction (1j]=1 and 12=0),equation 


(2.43) reduces to 


(Diik® + Ass — Iw? )(Desk*® + Ass—I@* )(Assk* —p w*) 
=) Drak) ( Assk sap w )—(Assk)\*(Desk>+Asas—im*)=0. (2.44) 
For waves propagating in the y-direction, equation (2.43) reduces to 
(Do2k* + Asa—Ia* )(Desk* + Ass —la*)(Aask® —p w°) 
—(Disk? )?(Assk® —p w”)—(Assk)* (Desk? + Ass—lm*)=0 (2.45) 
For unidirectional or cross-ply laminates, equation (2.43) can be further 
simplified (since D,6=0) for waves propagating in principal directions. For 
waves propagating in the x-direction (2.44) reduces to 
(Disk* + Ass —la* )(Assk* — p w*)—Ass* k* =0 
(2.46) 
For waves propagating in the y-direction (3.45) reduces to 
(Dz2k* + Aas—l@* (Aask® —p w*)— Aas? k® = 0 
(2.47) 
It can be recognized that the characteristic equation has more than one root. 
However, only one root approaches zero circular frequency as the wave number 
approaches zero. This is the root corresponding to the flexural branch of the 


frequency spectrum for plate waves. 
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lil. DISPERSION MEASUREMENTS 


A. INSTRUMENTATION AND EXPERIMENTAL PROCEDURE 

Figure 3.1 is a schematic of the instrumentation used for the phase velocity 
measurements. An arbitrary function generator (AFG), (LeCroy model 9100) 
was uSed to generate a 10 volt, peak-to-peak, gated nine cycle sine wave tone 
burst, shown in figure 3.2. The repetition rate of the tone burst was controlled 
by a 20 MHz pulse generator (Wavetek model 145). The tone burst from the 
AFG was amplified using a direct coupled amplifier (Krohn-Hite model DCA- 
50) over the frequency range of 0-500 kHz. The amplified signal was in turn 
input into a matching transformer (Khron-Hite model MT-55). The signal was 
used to drive a piezoelectric transducer which was used to transmit the wave into 
the plate. Several different transmitter and receiver transducer types were used 
to increase the size of the flexural mode. This issue is discussed in depth in 
section C of this chapter. Both exciter and receiver transducers were coupled to 
the plate using vacuum grease. The received signal was then amplified 60 dB 
using a model 1220A Physical Acoustics Corporation (PAC) preamplifier , in 
which the filter had been modified for broadband operation. Finally, the output 
signal was routed to a LeCroy 9400A digital oscilloscope to capture and digitize 
the waveform detected by the receiving transducer. 

In order to measure the phase velocity, a cursor on the oscilloscope screen 
was positioned at a reference point on the received wave. The receiving 
transducer was then moved a known distance, |. By maintaining the cursor at the 
reference point (phase point) of the waveform as the receiving transducer was 


translated, the time difference ,At, of the phase point was determined. Knowing 
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the time difference and the distance of transducer travel, it is a simple calculation 


to determine the phase velocity, Vp=I/At. 
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Figure 3.1 Schematic of Experimental Set-up 
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Figure 3.2 Input pulse 
B. EXPERIMENTAL RESULTS 

Shown in Figures 3.3 and 3.4 are the theoretical dispersion curves (from 
higher order plate theory) for the {O,90]4,5 laminate. Additionally, the 
experimental phase velocities are shown for the flexural waves propagating in 
the O and 90 degree directions. Figures 3.5 to 3.9 are the theoretical dispersion 
curves for the [0,45,90,-45]2,s laminate. The experimentally derived phase 
velocities are shown for the flexural waves propagating in the 0, 22.5, 45, 67.5 
and 90 degree direction. Material constants used for the calculation of the 
dispersion curves are given in Appendix A. Experimental data for all phase 
velocity measurements are given in Appendix B. 

The experimental data follows the general trend of theory. Except for the 
[O0,90]4 5 laminate measured in the 0 direction, the experimental phase velocities 
are consistently less than the theoretical. It is believed that this discrepancy 1s 


due to the actual elastic constants being less than those quoted by the 


manufacturer. 
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Figure 3.3 Theoretical Dispersion curves and data for [0,90 | 4.5 ,0 direction 
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Figure 3.4 Theoretical Dispersion curves and data for [0,90 |] 4.5 , 9% direction 
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Figure 3.5 Theoretical Dispersion curves and data for [0,45,90,-45 ] 2.5.0 direction 
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Figure 3.6 Theoretical Dispersion curves and data for [0,45,90,-45 | 2/5, 225° direction 
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Figure 3.7 Theoretical Dispersion curves and data for [0,45,90,-45 | 2.5 , 45 direction 
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Figure 3.8 Theoretical Dispersion curves and data for [0,45,90,-45 | 2.5, 67.5 direction 
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Figure 3.9 Theoretical Dispersion curves and data for [0,45,90,-45 ] 2.5 . 90 direction 
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Several subtleties associated with experimental phase velocity measurements 
were discovered which are discussed below. 

The selection of the transducer (size and resonant frequency) pulser/ receiver 
combination had an important result on the mode of the output wave. For 
example, when measuring the [0,90]4,s plate at the lower frequencies the 
combination of transducer A/B (as defined in Appendix B) produced an excellent 
presentation of the flexural mode (Figure 3.10). As the frequency is increased 
the extensional mode is introduced into the laminate (Figure 3.11). By switching 
the pulser/ receiver combination to B/C, the effect of the extensional mode is 
decreased (Figure 3.12). This may be due to the smaller diameter transducers 
not being excited by the extensional mode at this frequency. This is because the 
wavelength was much larger than the transducer diameter. The separation of the 
flexural and extensional mode ensured that the proper flexural mode was utilized 


to measure the phase velocity. 
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Figure 3.10 Reference wave form , (0,90]4,5 laminate, (frequency= 60kHz) 
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Figure 3.11 Extensional and flexural wave excitation, [0,90]4, laminate, 


(frequency= 130 kHz), A/B Transducer Combination 
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Figure 3.12 Delay of extensional mode due to transducer switch, [0,90]4,s 


laminate ,130 kHz, B/C transducer combination, 


Delaying the onset of the extensional mode relative to the flexural mode was 
also important to eliminate the problem of the extensional mode mixing and 
interfering with the flexural mode. This is a significant problem because if the 
two modes mix an incorrect phase velocity would be measured. 

The upper limit of flexural mode phase velocity measurements was realized 
when the extensional mode interfered and overpowered the flexural modes. For 
the [0,90]4,s plate this occurred at 220 kHz in the 90 direction and 260 kHz in 
the O direction. For the [0,45,90,-45]2,s laminate the upper limit was generally 
realized at 200 kHz for all directions. 

Because the experimental laminates were of finite lengths, wave reflections at 
certain frequencies were of concern. The mixing of the signal and reflections 
could produce an erroneous phase velocity. Because it was not possible to cancel 
all reflections, some frequencies were not measured due to this phenomenon. 

In summary, the experimental measurements of phase velocity required an 
understanding of the wave propagation. It was imperative to differentiate 
between the flexural and extensional modes and to realize the effects of the 
retlections. 

The first aspect of this thesis was concerned with verifying the higher order 
plate theory by comparing experimentally measured phase velocities to 
theoretically determined phase velocities. The second aspect of this thesis will 
now be concerned with using the experimental phase velocities to determine the 
laminate material properties by inverting the governing equations. This will be 


discussed in the following section. 
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IV. ALGORITHMS FOR RECOVERING MATERIAL 
CONSTANTS 


A. [0]14g LAMINATE 

Figure 4.1 is a comparison of dispersion curves using both classical plate 
theory (CPT) and higher order plate theory (HOPT) for the [O]j6 laminate. 
There is agreement, within a 5% variation, between the two theories if the 
frequency is limited to 19.0 kHz. At frequencies higher than 19.0 kHz, greater 
error arises in the phase velocities predicted by CPT due mainly to neglecting the 


shear deformation and rotary inertia. 
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Figure 4.1 CPT vs. HOPT [0],6 laminate , 0 direction. 


Rearranging equation (2.13) to solve for Dj; yields 


Vp*oh 
Diz 4.1 
1 Wore: ) ( ) 





Consequently, by using classical plate theory at the lower frequency regime, 
it is possible to recover Dj; by measuring the experimental phase velocity along 
the x-direction. The experimental frequency, plate thickness and plate density 
are known. Table 4.1 is a comparison of experimental values of D1, against the 
theoretical values. Classical lamination theory, which assumes static loading, was 
used to calculate D,,, All experimental numbers used in the calculations of 
material properties were described in the previous chapter. 

Different frequencies were investigated to determine the effect of frequency 
on the maternal properties. The experimental error at 6 kHz was greater due to 
higher uncertainties when measuring the time differences due to lower signal to 
noise ratios, excitation of the transducer's resonant frequency and reflections. 
These phenomena are discussed in mure detail at the end of this section. From 
analyzing the 8 kHz and the 10 kHz measurements, the values for D;; were found 
to be consistent within experimental error. Further analysis should include 
multiple experimental measurements at the same frequency to differentiate 
between systematic and random errors. 

Experimental errors were determined from procedures set forth by Holman 
[Ref. 15]. Appendix C shows all calculations and experimental uncertainties 
used for the determination of the total experimental error. 


Table 4.1 [O]16 laminate, Experimental and theoretical Dj 


Frequency (kHz) D)) (Nm) Di; (Nm) 
(Experimental) (Laminate Plate 
Theory) 
6.0 148.1 (4+/- 9.62) 140.0 
8.0 120.1 (+/- 8.26) 140.0 


10.0 ie C= S07) 140.0 

The discrepancy between experimental results and laminate plate theory is 
most likely due to the typical material constants supplied by the manufacturer 
being greater than the “as-manufactured" material constants of this particular 
plate. An area for further analysis would be to compare acoustically measured 
material properties to those derived from mechanical tests of specimens from the 
experimental laminates. 

Figure 4.2 is a graphical comparison of dispersion curves using CPT and 


HOPT of the [0]16 laminate for waves propagating in the 90 degree direction. 
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Figure 4.2 CPT vs. HOPT [O];6 laminate , 90 direction 


Rearranging equation (2.14) to solve for D22 yields 


Vp" ph 


Dye ead 
a PP 


). (4.2) 


when measuring phase velocities of flexural waves propagating in the y- 
direction. 

As before in the x-direction, by limiting experimental frequencies to the 
regions where classical plate theory and higher order plate theory agree, it was 
possible to use CPT to recover Do2. Table 4.2 is a comparison of experimental 


D2 against the theoretical D2 . 


Table 4.2 [O]16 laminate Experimental and Theoretical D2 


Frequency (kHz) D222 (Nm) D22 (Nm) 
(experimental) (plate theory) 

6.0 | 10.4 (4/- 1.38) o> 

8.0 P32 Ge 220) 935 

10.0 12.2 (+/- 1.97) 9:35 


The experimental results are within experimental error of each other but the 
8 and 10 kHz measurements were not within experimental error of laminate plate 
theory. Again, an area for further analysis would be to compare acoustically 
measured material properties to those derived from mechanical tests of 
specimens from the experimental laminate. 

All remaining material properties were calculated using higher order plate 
theory. 


Equation (2.46) algebraically rearranged to solve for Ass5 yields 
sae. Oo. Gyewe 
Ag, = 2 OS OCS 
Diik ep ORE, —lw*k 


Likewise, equation (2.47) algebraically rearranged to solve for A44 yields 
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ee p w*(D22k* —Ia*) 
D22k* —p w* —lw*k* 


(4.4) 

The first attempt at solving Ag4 and As5 was made using the lower 
frequencies and phase velocities used to calculate D}; and D2. This first attempt 
yielded negative values of Agg and As5 which was obviously an erroneous result. 
This error was due to the fact that (4.3) and (4.4) are quadratic equations with 
more than one root. The experimental frequency was increased to excite the 
positive root of Aqg4. Since this was only a quick check of an algorithm to 
recover material constants, only A4g4 was investigated for an experimental 
recovery. The frequency was increased to 100 kHz and the subsequent measured 
phase velocity resulted in an Agq of 5.38 MPa which compares favorably, within 
experimental error, to the theoretical value of 6.78 Mpa. 

Since Dy6 = O for a unidirectional plate, the following equation can be 
factored from equation (3.44): | 

Desk* + Asaa—Iw* = 0. (4.5) 


Equation (5.5) rearranged to solve for D66 gives 
Us= ley! — Ace (4.6) 

ke 
This is a quadratic equation in terms of the frequency. Recovering the 
positive root of this equation dictated the experimental frequency be increased to 
over 300 kHz. This presented a serious problem since in this frequency regime 
the extensional mode dominates, making the experimental measurement of the 
flexural mode impossible. Consequently, it was not possible to recover Dee and 
ultimately D 2. Dj 2 was to be calculated by equation (2.43) when measuring 


phase velocities in off-principal directions with all other material properties 


known. 


Sl 


In summary, by restricting the experimental frequency to the regime where 
CPT and HOPT agree, values of D,; and Do2 were recovered from CPT. 
However, measuring at these lower regimes presented their own problems and 
errors and was found to be difficult. The transducers used were not very 
sensitive at these lower frequencies resulting in low signal to noise ratios. The 
unidirectional laminate was especially prone to reflections when measuring in the 
90 degree direction. The waves propagated significantly faster in the 0 degree 
direction than in 90 degree direction. Consequently, when measuring in the 90 
direction, the reflections of waves propagating in the 0 direction interfered 
with the direct wave propagation along the 90 direction likewise making 
measurements difficult. The PAC R15 transducer was utilized for measurements 
in the O direction. When measurements were accomplished with the PAC R15, 
higher frequencies were introduced due to the input gated sine wave. The higher 
frequencies excited the transducer's resonant frequency which "rode" on top of 


the input frequency making measurements difficult. 


B. [0,45,90,-45]2.5 LAMINATE 

While employing classical plate theory for the recovery of Dj; and D292 for 
the unidirectional plate was successful, the inability to recover D2 and Dé6 
dictated the use of a more powerful algorithm. For the [0,45,90,-45]2 5 laminate, 
it was not possible to calculate Aq4g and As55 directly as before because D6 was no 
longer zero. It was determined that a desirable algorithm would be based upon 
equation (2.44). By assuming wave propagation along principal directions (0 or 
90 degrees) a significant simplification of the governing equation results. By 
successfully incorporating classical plate theory to determine Dj; and D9, the 


unknown material proj erties in (2.44) were reduced to 4 (Aga, As55, Die and 


oF 


D66). Once these unknowns were accounted for, it would then be possible to 
determine D2 from equation (2.43). 

CPT and HOPT agree within 5% if the experimental frequency is limited to 
16.5 kHz when measuring in the 0 direction (Figure 4.3). Equation (4.1) was 
applied to recover Dj; when limiting the experimental frequency to less than 
16.5 kHz. Table 4.3 1s presented as a listing of the experimentally derived values 


of Dj); against laminated plate theory. 
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Figure 4.3 CPT vs. HOPT [0,45,90,-45]2 ; laminate , 0 direction 
The measurement at 10 kHz was not within experimental error of plate 
theory. However, the remaining measurements agree with theory within 


experimental error. To determine if the 10 kHz measurement represents a 


os 


Systematic or random error, future studies should include statistical analyses 
using multiple measurements at the same frequency. 


Table 4.3 [ 0,45,90,-45]2; laminate, Experimental and Theoretical Dj, 


Frequency (kHz) Di, (Nm) Di, (Nm) 
(experimental) (plate theory) 
6.0 83.6 (+/- 9.63) 84.3 
0) 86.0 (+/- 6.43) 84.3 
8.0 ee Gqi= See) 84.3 
oa) 90.7 (+/- 7.14) 84.3 
OY, 66.8 (+/- 8.14) 84.3 


CPT and HOPT agree within 5% if the experimental frequency is limited to 
20.0 kHz when measuring in the 90 direction (Figure 4.4). Equation (4.2) was 
applied to recover D22 when limiting the experimental frequency to less than 
20.0 kHz. Table 4.4 is presented as a listing of the experimentally derived 
values of D22 against laminate plate theory. 

The experimental D2? was consistently lower than theory. As with the 
unidirectional laminate, this can be explained by the manufacturer's supplied 
material constants being higher than the as-manufactured experimental laminates 
material properties. The measurement at 5 kHz was not within experimental 
error of the other measurements. Multiple measurements at the same frequency 
should be made in future studies to distinguish between systematic and random 


Crro;rs. 
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Figure 4.4 CPT vs. HOPT, [0,45,90,-45]2. laminate , 90 direction 


Table 4.4 [ 0,45,90,-45]2; laminate, Experimental and Theoretical D22 


Frequency (kHz) D22 (Nm) D22 (Nm) 
(ex perimental) (plate theory) 
4.0 42.9 (+/- 5.02) S)e oe, 
5.0 51.88 (+/- 4.84) DS.> 
6.0 44.56 (+/- 4.77) 58.3 
7.0 41.0 (+/- 4.48) Der 
8.0 AS 2 Gi- 4.92) Des, 
9.0 40.5 (+/- 4.72) Sie, 
10.0 45.0 (+/- 5.84) Sieio3 


oe 


In order to solve for the four unknown material properties (A44, A55, D6 
and Dee), four measurements at four different frequencies in the O direction 
were made, thus resulting in four independent equations. Since equation 2.44 is 
nonlinear in terms of the material constants, a method to solve nonlinear 
equations was required. Wheatley [Ref. 16] describes how to reduce a nonlinear 
equation to a linear equation using a Newton's method approach which is 
summarized below. 

We begin our discussion of non-linear systems by assuming two functions 

t(x,y)=0, 
g(x,y)=0. 

Let x=r, y=s be a root. Both functions expanded as a Taylors series about 
the point (x;,y;) in terms of (r-xi),(S-yi), where (xi,y1) is a point near the root, 
gives 

f(r,s)=O= f(xi,y)+ fx(xi,yi)G-xi)+ fy, yi)(s-yi)t ... 
g(r,s)=O= g(Xj,yi)t+ Bx(Xi,yi(T-xXi)+Zy(Xi,yi(S-yI)+ ... (4.7) 
The subscript notation designates the partial derivatives. Truncating the 
series after the linear terms gives 
: pee | ‘| fe(xi, Vi) ty (Xi, yi) i Ml! (4.8) 
O} LOX Y)] LOY) GX, Vi) ILS— yi 
The above equation can be rewritten to solve as the system of equations 


~~. oda | ie (4.9) 
Qx( Xi, Yi) Qy(Xi, Yi) || Ayi g(Xi, yi) 


where Axi=r-xi and Ay; =s-yl. 
Equation (4.9) is solved by Gaussian elimination and an improved estimate of 


the root 1s made by setting 


dep (4.10) 
Yir+i yi AYi 
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The process is then iterated upon unul f and g are close to zero. 

A computer program was written in the BASIC language based upon the 
methods described by Wheatley. 

However, this method never converged towards a solution. The non- 
convergence was due to the ill-conditioning of the partial derivative, square 
matrix of equation (5.9) as a result of the large magnitude differences between 
Aa4, Ass (9.42 x 10° Pa m) and Dj6,D66 ( 19.9 Nm). The values of Aga4/A55 


were of a magnitude of a million times greater than the values of D16/Dée. 


C. [0,90]45 LAMINATE 

It was initially believed that for the cross-ply laminate at the experimental 
frequency required to recover Dee, the flexural mode would not be overpowered 
by the extensional mode. This was the limitation for the unidirectional laminate. 
The cross-ply laminate was analyzed next using the algorithm developed for the 
unidirectional plate since D,¢6=0. That is, using CPT in the lower frequency 
regimes where CPT and HOPT agree to recover Dj; and D272. HOPT is then 
used to recover the additional material constants. 

CPT and HOPT agree within 5% if the experimental frequency is limited to 
18 kHz when measuring in the 0 direction (Figure 4.5). 

Equation (4.1) was applied to recover Dj; when limiting the experimental 
frequency to less than 18 kHz. Table 4.5 is presented as a listing of the 


experimentally derived values of Dj; against laminate plate theory. 
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Figure 4.5 CPT vs. HOPT [0,90]4, laminate , 0 direction 


The measurements 8 and 14 kHz were not within experimental error of 
measurements made at the other frequencies. Statistical analyses using multiple 
measurements 1S necessary to distinguish between systematic and random errors. 


Table 4.5 [ 0,90]4., laminate, Experimental and Theoretical Dj 


Frequency (kHz) Di, (Nm) Di, (Nm) 
(experimental) (plate theory) 
8.0 60.3 (+/- 6.73) PSO 
10.0 92.3 (4+/- 9.66) 113.0 
12.0 88.6 (+/- 11.97) L130 


38 


14.0 69.9 (+/- 9.5) ela) 
16.0 83.3 (+/- 12.34) 113.0 
CPT and HOPT agree within 5% if the experimental frequency is limited to 


16.5 kHz when measuring in the 90 degree direction (Figure 4.6) 
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Figure 4.6 CPT vs. HOPT , [0,90]4.; laminate , 90 direction 


Equation (4.2) was applied to recover D22 by limiting the experimental 
frequency to less than 16.5 kHz. Table 4.6 is presented as a listing of the 
experimentally derived values of D2 against laminate plate theory. 


Table 4.6 [ 0,90]4, laminate, Experimental and Theoretical D2 


Frequency (kHz) D2 (Nm) D22 (Nm) 
(experimental) (plate theory) 
8.0 80.9 (+/- 9.45) 0) ao 
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10.0 51.2 (4/- 5.98) $1.3 


126 73.4 (+/- 9.60) S13 
14.0 GOG7aGe/- 9-05) 81.3 
16.0 $3.3 G/- 12.34) Slee) 


The measurement at 10 kHz was not within experimental error of the other 
measurements. Again, Statistical analyses is necessary to differentiate between 
systematic and random errors. 

From HOPT, equations (4.3) and (4.4) were applied to solve for As55 and 
A44. The average experimental A55 was calculated to be 9.83 Mpa-m which 
compares nearly exactly with the theoretical value of 9.85 MPa-m. The average 
experimental A4gq4 was calculated to be 7.80 MPa-m which is 20.8 % lower than 
the theoretical value of 9.85 MPa-m. Table 4.7 is presented as a listing of the 
experimentally derived values of As5s5 against theory while Table 4.8 is the 
experimentally derived values of Ag4 against theory. The small variation of the 
experimental Aqq and As55 can be attributed to the less difficult measurement 
techniques required when measuring at the higher frequencies. 

Equation (5.6) was utilized to attempt the recovery of D66. However, 
similarly to the unidirectional laminate, when the experimental frequency was 
increased to excite the positive root of (5.6), the extensional mode dominated, 
making the measurement of the flexural mode impossible. Subsequently, Dée 


and ultimately Dj2 were not recoverable from this procedure. 
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Frequency (kHz) 


60.0 

80.0 

100.0 
120.0 
130.0 
140.0 
150.0 
160.0 
170.0 
180.0 
1oU0 
200.0 
220.0 
240.0 
260.0 


Ass (MPa-m) 
(Experimental) 
9.82 (+/- 1.07) 
8.17 (+/-1.65) 
7.78 (t/-) 1,96) 
9.89 (+/- 1.24) 
10.11 (+/- 1.25) 
10.33 (+/- 1.26) 
eG) = | IG) 
10.47 (+/- 1.23) 
9.94 (4/- 1.11) 
Zio) A Sel ead wa Uy) 
7 Ge lew) 
10.04 (+/- 1.17) 
POG) eas) 
ioe GeO 1) 
9.84 (+/- 1.67) 


Table 4.7 [ 0,90]4.; laminate, Experimental and Theoretical A55 


Ass (MPa-m) 
(plate theory) 
9.85 
9.85 
DES 
9.85 
9.85 
Des 
N85 
Shots) 

DES 
285 
Dees 
si 
9.85 
9.85 
ISS 


Table 4.8 [ 0,90]4, laminate, Experimental and Theoretical Agq4 


Frequency (kHz) 


60.0 
80.0 
100.0 


Aaa (MPa-m) 
(Experimental) 
Scolar ESI): 
7.26 (+/- 0.916) 
ra U2) 


4] 


A44 (MPa-m) 
(plate theory) 
5/10) 

9.85 
9.85 


110.0 7.39 (+/- 0.839) 9.85 
120.0 7.66 (+/- 0.783) 9.85 
130.0 7.74 (+/- 0.854) 9.85 
140.0 7.97 (+/- 0.874) 9.85 
150.0 7.86 (+/- 0.699) 9.85 
160.0 8.01 (+/- 0.858) 9.85 
170.0 8.02 (+/- 0.808) 9.85 
180.0 8.00 (+/- 1.34) 9.85 
190.0 TG (EE DS 9.85 
200.0 7.90 (+/- 1.58) 9.85 
210.0 7.69 (+/- 1.51) 9.85 
220.0 7.8 (+/- 0.97) 9.85 


D. DISCUSSION / SUMMARY OF ALGORITHMS TO RECOVER 
MATERIAL CONSTANTS 


When first reviewing the procedures to recover the material constants, it was 
initially believed that the algorithm used for the unidirectional and cross-ply 
laminates would be sufficient. The employment of classical plate theory to 
recover D,,; and D2? was attempted in all three laminates with limited success. 
Further work should include multiple measurements at the same frequency so 
statistical analyses can distinguish between systematic and random errors. The 
utilization of HOPT was successful to recover the transverse extensional terms 
A44 and As55 for the unidirectional and cross-ply laminates. The experimentally 
derived material properties were generally lower than values predicted by 
laminate plate theory. This is consistent with the experimentally measured phase 


velocities being generally lower than phase velocities predicted by HOPT. 
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Further analyses should also include comparisons between acoustically measured 
material properties to those derived from mechanical tests of specimens from the 
experimental laminates. The recovery of Dee dictated experimental frequencies 
where the extensional mode dominated making the measurement of the flexural 
mode impossible. 

If Deé6 could have been recovered, the only material constant not accounted 
for would have been Dj2. It would have possible to solve for Dj2 by the 
substitution of all known material constants and an off principal axis velocity 
measurements into equation (3.43). 

The inability to recover Dg66 and Dj2 in the cross-ply and unidirectional 
laminates dictated the exploration of a new procedure to recover the material 
constants in the angle-ply laminate. It was determined that a new algorithm 
would be based upon equation (2.44), (2.44) was chosen because a significant 
simplification of the governing equation ensues when measuring wave 
propagation in the principle directions (Oor 90). By successfully 
incorporating CPT at the lower frequency regimes, it was possible to recover 
Dj; and D22. The unknown material properties in (3.44) were reduced to 4 
(Aga, Ass, Di6 and Dee). Aga and Ass can not be calculated directly as before 
because Dj 6 is not zero. 

In order to solve for the four unknowns, four experimental conditions were 
input to produce four independent equations. These independent equations are 
nonlinear in terms of the material constants. A numerical method developed by 
Wheatley [Ref. 15] was then utilized to reduce a nonlinear equation to a linear 
equation using a Newton's method approach. Unfortunately, this method never 


converged towards a solution. The non-convergence was due to the ill- 
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conditioning of the partial derivative, square matrix of equation (4.9). It was 
later discovered that the large magnitude differences between A4g4/ A55 and 
D16/D66 is responsible for the ill-conditioning of the matrix. 

Although this method base upon Wheatley did not converge, it 1s believed 
that a solution can be derived from a more robust numerical method. When 
factoring out equation (2.44) or (2.45), it can be proven that the equation is 
nearly quadratic in terms of the material properties. Reference 16 recommends 
several numerical methods to undertake to solve more advanced problems of this 
type. A methods which show promise is the Conjugate Direction method for non- 
linear systems. 

Another method to attempt to solve for the material constants which shows 
promise but was not attempted due to time constraints is based upon equation 
(2.9) using classical plate theory. (2.9) can be solved using linear algebra since 
the equations are linear in terms of the material constants. By measuring phase 
velocities along principal directions in the lower frequency regimes, we have 
proven it is possible to recover D,,; and Do». The unknown material properties 
of equation (2.9) are reduced to three ,D12, D16 =D26 (symmetric laminate only) 
and De66. By measuring in off-principal directions at varying experimental 
frequencies, it will be possible to generate three linearly independent equations in 
terms of the unknown properties. This is the very approach that was attempted 
in the [0, 45,90,-45]4.; laminate that did not converge due to ill-conditioning. 
However, this situation is different in that all unknown parameters are roughly 
of the same magnitude which would prevent the ill-conditioning that occurred 
before. Also this numerical method represents a system of linear equations 


which is a significant simplification over a system of nonlinear systems. 
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V. SUMMARY 


Composite materials promise to be a major structural material of the future. 
Composite materials present unique problems in that the complexities of the 
manufacturing process present increased opportunities for the introduction of 
defects. An ability to nondestructively determine the elastic properties of the 
composite component would be an invaluable tool for quality assurance. This 
thesis investigated a technique to yield the stiffness constants of composite 
laminates. 

The initial aspect of this thesis was the verification of flexural wave 
propagation theory in composite laminates. Experimental phase velocities were 
measured in graphite/epoxy plates and were compared with dispersion curves 
calculated using higher order plate theory. Figures 3.3 through 3.9 summarize 
the results. The general trend of theory is observed but the phase velocities are 
typically lower than those predicted theoretically. This discrepancy could be 
explained if the actual elastic constants of these plates were smaller than those 
quoted by the manufacturer. An area for future analysis would be to compare 
acoustically measured material properties to those derived from mechanical tests 
of specimens from the experimental laminates. 

The second aspect of this thesis was concerned with using the experimental 
phase velocities to determine the laminate material properties by inverting the 
governing equations. This turned out to be a very difficult problem for several 
reasons and different methods were tried. The first method incorporated 
Classical Plate Theory (CPT) to recover D11 and D22 and Higher Order Plate 


Theory (HOPT) to recover A44 and A55 for the unidirectional and cross-ply 
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laminates. However, at the higher frequencies necessary to recover D66 and D12 
the extensional mode dominated making the measurement of the, flexural mode 
impossible. 

The inability to recover D66 and D12 dictated the exploration of a new 
procedure on the [0,45,90,-45]2,s laminate. CPT was applied to recover Di1 and 
D22. By assuming wave propagation along a principal direction, the wave 
propagation equations were significantly simplified. The unknown material 
properties were reduced to four. Four experimental velocity measurements at 
different frequencies along the zero degree direction were input to produce four 
independent equations. The four equations can be solved for the four remaining 
unknown parameters. Because these equations are nonlinear in terms of the 
material constants, a nonlinear numerical method was utilized. Unfortunately, 
this method did not converge towards a solution due to the ill-conditioning of the 
matrix. Several other numerical methods were suggested in Chapter IV for 
further exploration. 

A method proposed but not examined was to use non-principal direction 
velocity measurements in conjunction with classical plate theory (developed in 
Chapter 2) to recover D12, D16 and D66. These measurements must be made at 
frequencies at which classical plate theory is valid. Previous measurements at the 
lower frequency regime have proven to be difficult but the numerical method to 
recover the unknown parameters would be greatly simplified because the 


equations are linear. 
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APPENDIX A- LAMINATE MATERIAL PROPERTIES 


The dispersion curves were calculated using the following material 
properties of AS/4 3502: 

E)1;=144.8 GPa  E2=9.65 GPa G)2=5.97 GPa 

Vi2=0.3 V23=0.34 . 

All plates were assumed to have a density of 1550 kg/m?. 

The thickness of the [O]16 plate was 0.002352 m. 

The thickness of the [0,90]4 5 plate was 0.0024688 m. 

The thickness of the [0,45,90,-45]2., plate was 0.0023368 m. 

The calculated bending stiffnesses, Diy, for the [0]16 laminate are : 

i= 115-0 Nm, Dy9=3.65 Nm, Dig= Do6= 0.0 Nm, 

D22=81.4 Nm, Dee=7.49 Nm. 

The calculated transverse extensional stiffnesses, Aij, for the [0,90]4., 
laminate are A44=9.84Mpa. and A55=9.84 MPa. 

The calculated bending stiffnesses, Diy, for the [0,90]q.; laminate are 

D,,=113.0 Nm ,D)2=3.65 Nm, D16= D26= 0.0 Nm, 

D79=81.4 Nm, Deé6=7.49 Nm. 

The calculated transverse extensional stiffnesses, Aj, for the [0,90]4.; 
laminate are A44=9.84Mpa. and A55=9.84 MPa. 

The calculated bending stiffnesses , Diy, for the (0,45,90,-45]2., laminate are : 

D11=84.3 Nm D)2=1.67Nm Dj¢6= Dog = 5.93 Nm 

D22=53.8 Nm Deée=19.9 Nm. 

The calculated transverse extensional stiffnesses , Aij, for the [0,45,90,- 


45 ]25 laminate are A44=9.42MPa and A55= 9.42 MPa. 


4] 


The shear correction factors for composite laminates, kiki were 


experimentally found to be 5/6 by Tang, Stuffler and Henneke [Ref. 12]. 
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APPENDIX B - EXPERIMENTAL DATA 


The transducers used for experimental measurements were: 


A- Harisonic HC-483, resonant frequency 2.25 MHz, diameter 12.7 mm 
B- Harisonic HC-483, resonant frequency 1.0 MHz, diameter 12.7 mm 


C- Harisonic GO504 , resonant frequency 5.0 MHz, diameter 6.3 mm 


D- Physical Acoustics Corp R15-3194, resonant frequency, diameter 12.7 


E- Accelerometer, diameter 12.7 mm. 
N.A.- Not Available 


Table 1. Experimental Phase Velocity Measurements 


[0,90]4.5 laminate , 0 propagation 


| Frequency 


(kHz) 


8.0 
10.0 
12.0 
14.0 
16.0 
60.0 
80.0 

100.0 
120.0 


distance, 


(m) 


0:1 
0.1 
0.1 
0.13 
0.1 
Onl 
0.1 
0.1 
0.1 
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time 


difference 


(US) 
224 
180 
166 


Phase 


Velocity 


(m/s) 
446.4 
SDR: 
602.4 
Gle.2 
684.9 
1116.1 
1179.2 
1314.0 
| as 


oie) 
Receiver 
Transducer 

A/B 
A/B 
A/B 
A/B 
A/B 
A/B 
A/B 
A/B 
A/B 


SO 


130.0 0.1 iO 1408.4 C/C 
140.0 0.1 69.4 1440.9 B/C 
150.0 0.1 69.2 1445.0 B/C 
160.0 0.1 67.2 1488.1 B/C 
170.0 0.1 67.6 1479.3 C/C 
180.0 0.1 68.4 1462.0 C/C 
190.0 0.1 66.4 1506.0 C/C 
200.0 0.05 32.8 1524.0 C/C 
220.0 0.06 39.4 S22 ae C/C 
240.0 0.1 64.6 1548.0 Ge 
260.0 0.06 38.4 1562.5 One 
Table 2. Experimental Phase Velocity Measurements 
[0,90]4.; laminate ,90 propagation 
Frequency distance, | time Phase Exciteg, 
(kHz) (m) difference Velocity Receiver 
(us) (m/s) Transducer 
8.0 0.05 104.0 480.7 A/B 
10.0 0.07 146.0 479.5 A/B 
12.0 7 174.0 574.7 A/B 
14.0 0.1 165.0 606.0 A/B 
16.0 IGS 720 684.9 A/B 
60.0 0.1 93.6 1068.4 A/B 


80.0 0.1 88.4 Sl 2 A/B 


100.0 0.1 Boo 1201.9 A/B 
110.0 Or $1.6 25> A/B 
120.0 Cala Nye k2 61.5 A/B 
130.0 Or 77.8 S555 A/B 
140.0 0.1 76.0 115.6 A/B 
150.0 wale 90.6 1324.5 A/B 
160.0 0.1 74.0 351.3 A/B 
170.0 0.1 74.2 1347.7 Ge 
180.0 0.06 44.1 ISO Cie 
190.0 0.06 43.8 1369.8 Ce 
200.0 0.05 36.2 i339 C/C 
210.0 0.05 36.4 [37376 Cie 
220.0 0.08 57.4 SSH Ce 


Table 3 Experimental Phase Velocity Measurements 


[0,45,90,-45]4, laminate , 0 propagation 


Frequency distance, | time Phase Ee XCiten, 
(kHz) (m) difference Velocity Receiver 
(us) (m/s) Transducer 
20.0 0.1 157.0 700.0 N.A 
30.0 0.1 151.0 $27.0 N.A 
40.0 al 103.0 970.0 N.A 
60.0 0.1 cles 1096.0 N.A 


al 


80.0 0.1 82.8 1207.8 
100.0 0.1 78.0 128200 
120.0 0.1 76.0 [31S 
140.0 0.1 Tae? 1347.7 
160.0 0.09 65.0 1384.6 
180.0 0.1 70.0 1428.0 
200.0 0.1 50.0 1500.0 


Table 4 Experimental Phase Velocity Measurements 


(0,45,90,-45]4.; laminate , 22.5 propagation 


Frequency distance] —st time Phase 
(kHz) (m) difference Velocity 
(LS ) (m/s) 
40.0 0.1 3.0 884.9 
60.0 0.1 98.0 1020.4 
80.0 0.1 85.6 1168.2 
100.0 0.1 $2.0 I ZAMEL SS 
120.0 0.1 76.6 1305.4 
140.0 0.1 WOs0 Noe 58 
160.0 0.1 G22 1385.0 
180.0 0.08 56.6 1413.4 
200.0 0.08 54.8 1459.8 


N.A 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 


Excitem, 
Receiver 
Transducer 
N.A 

N.A 

N.A 

N.A 

N.A 

N.A 

N.A 

N.A 

N.A 


Table 5 Experimental Phase Velocity Measurements 


[0,45,90,-45]4.; laminate , 45.0 propagation 


Frequency 


(kHz) 


30.0 
40.0 
60.0 
80.0 
100.0 
120.0 
140.0 
160.0 
180.0 


distance, | 


0.1 
0.1 
0.1 
0.1 


0.1 


0.1 
0.1 
0.1 
0.1 


time 


difference 


( [LS ) 
128.0 
119.0 
100.8 
91.2 
84.8 
SD 
77.8 
74.8 
Fa 


Phase 
Velocity 
(m/s) 
TASS: 
840.3 
eal 
1096.5 
LS 
L234 
35.5 
| By SiOn 
1366.1 


Table 6 Experimental Phase Velocity Measurements 


(0,45,90,-45]4.; laminate , 67.5 propagation 


Frequency 


(kHz) 


20.0 
40.0 


distance, | 


0.1 
0.1 


time 


ditference 


(US ) 
160.0 
120.0 


Phase 
Velocity 
(m/s) 
625.0 
So Iee, 


E Citely 
INcCCelt el 
Transducer 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 


BxCiicy, 
Receiver 
Transducer 
N.A 

N.A 


50.0 
60.0 
80.0 
100.0 
120.0 
140.0 
160.0 
180.0 


0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
On 
0.1 


110.8 
102.8 
96.0 
86.6 
$1.8 
80.4 
78.6 
158 


SOV 

Zen 

1041.0 
Elo2al 
222 
1243.8 
1272.3 
IS Sige 


Table 7 Experimental Phase Velocity Measurements 


[0,45,90,-45]4., laminate, 90 propagation 


Frequency distance, | time Phase 
(kHz) (m) difference Velocity 
(US) (m/s) 
20.0 Onl 156.0 641.0 
30.0 0.1 [350 740.7 
40.0 0.1 ee) 826.4 
60.0 0.1 102.8 7 38 
80.0 0.1 92.4 IOSZe2 
100.0 0.1 $6.4 1157.4 
120.0 0.1 $2.8 WADE # 
140.0 0.1 78.2 Zio 
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N.A 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 
N.A 


Excite, 
Receiver 
Transducer 
N.A 
N.A 

N.A 

N.A 
N.A 

N.A 

N.A 

N.A 


160.0 0.1 TES 1288.0 
170.0 0.1 77.0 LAD 
180.0 0.1 12 78 


Table 8 Experimental Phase Velocity Measurements 


(O]i¢ laminate , O propagation 


Frequency distance, | time Phase 
(kHz) (m) difference Velocity 
(US) (m/s) 
6.0 0.2 404.0 495.0 
8.0 OZ 368.7 542.0 
10.0 0.2 330.0 606.0 


Table 9 Experimental Phase Velocity Measurements 


(O]i6 laminate , 90 propagation 


ned 
(KHz) 


6.0 
8.0 
10.0 
100.0 


distance, | 


(m) 


Onl 
0.1 
0.1 
0.1 


time 


difference - 


(US ) 
392.0 
320.0 
292.0 
113.6 


35 


Phase 


Velocity 


(m/s) 
Pos 
SAS 
342.5 
880.3 


N.A 
N.A 
N.A 


Exciter) 
Receiver 
Transducer 
B/D 

B/D 

B/D 


Exeiicn, 
Receiver 
Transducer 
B/E 

B/E 

B/E 

B/E 


APPENDIX C - EXPERIMENTAL ERRORS 


Let R be a function of some independently measured quantities xl, x2, 
X35. xe litatels 
R=R(x],X2,X30eee ne (C-1) 
Let Wry be the uncertainty in the result, R, and W1, W2, ... Wn be the 
uncertainties in the independent variables. Then the uncertainty in the result is 


given in Ref. [15] as 





/2 
Wa =| (S20)? + e)s eye Wo) | (C-2) 
OX: OX 2 OXn 
A. D11/ D22 
For this application, the governing equation was 
Vo"ph 
Du = (355) | (C-3) 


The phase velocity, Vp. was divided into the experimental measurements of 
distance and time difference , At. 
Adist 
Vo 
P At 





Ces) 


The density, p , was likewise divided into the experimental measurements of 


mass and volume (length, width and height). 
=o 2 ee 
volume = Iila2h 








Substituting experimental measurements (C-4) and (C-5) into (C-3) gives 
_Ve'ph ___(dist)*(kg)h 





DA= = . (C-6 
ane (At)* (42°F? )inlah .°) 
Canceling the height, h, results in 
Bees Oe (C-7) 


(At) (42°F? hile’ 
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The following experimental uncertainties were assumed when measuring at 
the lower frequency regimes: 

(A) dist=dist +/- 0.002m 

(B) kg= kg +/- 0.0001 kg 

me) At= At +/- 5yus 

(D) ly=l, +/- 0.001m 

(E) 1l2=12 +/- 0.001m 

(F) experimental error from input frequency was small in comparison to 
other experimental parameters and thus was neglected. 

After partial differentiation in terms of the independent variables, the 


uncertainty in the experimental result of D,; was 
Algisnieng (dist)° ; 
ee Aisne) 
A ae Phils At 4eP hls gy + 
(dist)" kg : (dist)’ kg 2 41/2 
= Al1)° + (—— Al (COs 
AC ZAG am ae i ip At ae hile ame” Sue) 


4(dist)* kg 


ieee Nt | 
NEA eof hile 


+( 


Substituting experimental data and uncertainties into (C-8) gives 
experimental errors for D11 quoted in Chapter IV. 
The same procedures and uncertainties were used for the calculation of 


experimental error in D722. 


B. A44/ A55 


The governing equation for the calculation of A55 was 
en op w* (Disk? —!@*) 


—— . C-9 
Duik* —p w* —lw*k? \ 


p was divided into the experimental measurements of mass (kg.) and the 


length and width of the plate resulting in 
Papp (C-10) 
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The inertia term, I, was similarly divided into the experimental 


measurements of mass, height, length and width, resulting in 


2 
ee (C-11) 
12/ 1/2 


(C-10), (C-11) and (C-4) were substituted into (C-9), resulting in the 





following experimental expression of AS55: 
16Diikgn*f*At*? —16kg*h* nt? 


22 ee : 
Ass = = latte OTD 
16D" At? 4kgn*f® —16Duikgh*n*f* at’ 
dist* Ih2 12izdist® 


In accordance with (C-2), (C-12) was differentiated with respect to the 
experimental parameters. All experimental uncertainties from D11 were 
maintained except that the time difference uncertainty, At, was reduced to de 
lus because the output signal was much clearer when measuring at the higher 
frequencies and the laminate height measurement uncertainty was small and thus 
neglected. Experimental data and uncertainties were input to determine 
experimental uncertainties quoted in Chapter IV. 

The same uncertainties and procedures were used to calculate the 


experimental uncertainties of A44. 
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